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ABSTRACT 



BL Lac objects are the best candidates to study the jet properties since their 
spectral energy distributions (SEDs) are less contaminated by the emission from 
the accretion disk and external Compton processes. We compile the broadband 
SEDs observed with Fermi/LAT and other instruments from literature for 24 TeV 
BL Lac objects. Two SEDs, which are identified as a low or high state according 
to its flux density at 1 TeV, are available for each of ten objects. These SEDs 
""oj can be explained well with the synchrotron+synchrotron-self-Compton model. 

We constrain the magnetic filed strength {B) and the Doppler factor (6) of the 
radiation region by incorporating the ^^-minimization technique and the 7-ray 
transparency condition. Twenty-four out of the 34 SEDs can constrain B and 6 
in 1 0" significance level and others only present the limits for both B and S. The 
typical values of B and 6 are 0.1 ~ 0.6 G and 10 ~ 35, respectively. The derived 
values of 7b are significantly different among sources and even among the low and 
high states of a given source. Prominent flux variations with a clear spectral shift 
are observed and a tentative correlation between the ratio of the flux density at 1 
TeV and the ratio of 7b in the low and high states is presented, likely indicating 
that the relativistic shocks for the electron acceleration may be responsible for 
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the flux variations and the spectral shift. A weak anti-correlation between the jet 
power and the mass of the central black hole is observed, i.e., Pjet oc M^^, which 
disfavors the scenario of a pure accretion-driven jet. Implications for the blazar 
sequence and the intergalactic magnetic field from our results are also briefly 
discussed. 

Subject headings: radiation mechanisms: non-thermal — BL Lacertae objects: 
general — gamma-rays: observations — gamma-rays: theory 

1. Introduction 

The broadband spectral energy distributions (SEDs) of blazars are double-peaked. Gen- 
erally, the bump at the IR-optical-UV band is explained with the synchrotron process of 
relativistic electrons accelerated in the jets, and the bump at the GeV-TeV gamma-ray band 
is due to the inverse Compton (IC) scattering of the same electron population (e.g., Ghis- 
ellini et al. 1996; Urry 1999). The seed photon field may be from the synchrotron radiation 
themselves (the so-called SSC model; Maraschi ct al. 1992, Ghisellini et al. 1996) or from 
external radiation fields (EC), such as the broad line region (BLR; Sikora et al. 1994), ac- 
cretion disk (Dernier ct al. 1992), torus (BLazejowski et al. 2000), and cosmic microwave 
background (CMB; Tavccchio et al. 2000). The low energy bump is well sampled with 
the data from multi-wavelength campaigns at the radio, optical, and X-ray bands. The 
observations with CGRO/EGRET and ground-based TeV gamma-ray telescopes sketch the 
feature of the high energy bump for some sources, but they are still poorly constrained. The 
observations with Fermi/LAT, which covers an energy band from 20 MeV to ~ 300 GeV, 
together with the ground-based observations at the TeV gamma-ray band, now pin down 
both the luminosity and the frequency of high energy peak well in the observed SEDs for 
some blazars. Broadband SEDs obtained simultaneously or quasi-simultaneously are critical 
to reveal the radiation mechanisms. 

About forty active galactic nuclei (AGNs) have been detected in the TeV gamma-ray 
band since the first detection of TeV gamma-rays from Mkn 421 with the Whipple imaging 
atmospheric-Cherenkov telescopes in 1992 (Punch et al. 1992) and most of them arc blazars. 
Observations show that these sources have violent variability in multi-frequency, especially 
at the X-ray and the TeV bands (Buckley et al. 1996; Takahashi et al. 2000; Sambruna et al. 
2000). The flux variations at the X-ray and TeV bands are correlated, such as that observed 
in Mkn 421 (Blazejowski et al. 2005). Some sources detected by EGRET are not confirmed 
by Fermi/LAT, and some in the IFGL (1st LAT Catalog) with fluxes much over the EGRET 
threshold were not detected by EGRET (Abdo et al. 2010a). Three TeV AGNs still have 
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not be detected by Fermi/LAT in the IFGL and in the 2FGL (Ackermann et al. 2011). 

These facts indicate that AGNs should have significant variabihty in the GeV-TeV band. In 
the high state of TeV emission, the peak frequency {ug) of the synchrotron radiation in the 
SEDs moves to higher energies for some sources, showing a tendency that a higher TeV flux 
corresponds to a harder spectrum for the emission at the X-ray and gamma-ray bands, as 
observed in Mkn 501 (Anderhub et al. 2009a). It is unclear what may be responsible for 
the flux variation and the corresponding spectral shift. The abundant data observed with 
Fermi/LAT in the MeV-GeV band, together with the TeV gamma-ray data, now provide an 
excellent opportunity to reveal the physical origin of the temporal and spectral variations of 
these TeV sources. 

It is generally believed that the jets of AGNs are powered by the central massive black 
hole (BH), whereas the mechanism of the jet production is still a mystery. The proposed 
models can be simplifled as two types, i.e., accretion and rotation of the BH. The connection 
between jet and accretion disk has been widely investigated (Ulvestad & Ho 2001; Marscher 
et al. 2002; Merloni et al. 2003; Falcke et al. 2004; King et al. 2011). With a sample of 23 
blazars in the three months of Fermi satellite survey, Ghisellini et al. (2009) found that the 
jet power is slightly larger than the disk luminosity and proportional to the mass accretion 
rate. However, the accretion process cannot explain the radio loudness of AGN, and thus 
the rotation of the BH is suggested to be the possible ingredient (Sikora et al. 2007; Lagos 
et al. 2009). In fact, there is some evidence for rapidly spinning BH (Cao & Li 2008; Wu 
et al. 2011). The jets of AGNs may be powered by both the accretion process and the 
spin of the central BH (Fanidakis et al. 2011). The SEDs of BL Lac objects dominated by 
jet emission suffer less contamination of the emission from accretion disk and EC process. 
Therefore, they are the best candidates for investigating the jet properties and the relation 
to the central BHs. 

In this paper, we compile a sample of SEDs for TeV BL Lac objects to study their 
physical properties of the jets and explore the physical reasons that result in the flux variation 
and the spectral shift in these sources as well as the relations of jet production to the accretion 
and rotation of the central BHs. The sample selection and the observed SEDs are presented in 
§2. The model, the strategy, and the results of SED fltting are described in §3. The physical 
properties derived from our SED fitting are presented in §4. A summary for our results is 
in §5. A simple discussion on implication of the blazar sequence and on the intergalactic 
magnetic field is given in §6. 
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2. Sample Selection 

All TeV BL Lac objects that have Fermi/LAT detections are included in our sample, 
except for six sources that have only upper limits of Fermi/LAT detections. Twenty-four 
BL Lac objects are included in our sample as listed in Table 1. We compile their broadband 
SEDs that were simultaneously or quasi-simultaneously observed with Fermi/LAT and other 
instruments from literature. Only one broadband SED is constructed for 14 objects and two 
broadband SEDs are available for the other ten objects. The two well-sampled SEDs for the 
ten sources are defined as a high or low state with the observed or extrapolated flux density 
at 1 TeV. The observed broadband SEDs are shown in Fig. [T] The description of each source 
is summarized in Appendix. 



3. Modeling the SEDs 
3.1. Model 

Although the contributions of EC component may be involved for some objects (e.g., 
3C 66A, Abdo et al. 2011), the external photon fields are usually much weaker than the 
synchrotron radiation photon field for the BL Lac object^. Multi-zone syn+SSC models are 
also used to explain the SEDs of some sources, such as Mkn 421 (Blazejowski et al. 2005). 
For seeking the consistency among objects in our statistical analysis, we only consider the 
simple single-zone syn-|-SSC model to fit the observed broadband SEDs in this paper. 

In the process of the SED modeling, the radiation region is taken as a homogeneous 
sphere with radius R. The electron distribution is assumed as a broken power law with 
indices pi and p2 below and above the break energy 7bmeC^, 



iV(7) = No 



7 7min < 7 < 7b, 

< ^ < 



where pi^2 = 2ai_2 + 1, and ai_2 are the observed spectral indices below and above the 
synchrotron radiation peak. The synchrotron radiation peak frequency of a single electron 
with Lorentz factor 7 is 

z/syn = 3.7 X lO^B^^-—. (2) 

i. -\- z 



^The fitting results with the SSC+EC model for some sources reported in the literatures are summarized 
in Appendix. 
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The synchrotron emissivity es(i^) is calculated with 

e,H = -/ djN{j)P,{u,j), (3) 

"'7min 

where Ps(z/, 7) is the single electron synchrotron emissivity averaged over an isotropic dis- 
tribution of pitch angles (see e.g., Crusius & Schlickeiser 1986; Ghisellini et al. 1988). The 
synchrotron radiation field /s(i^) is calculated by the transfer equation, 

W = '-^[l - e-'^'^n (4) 

where k{i') is the absorption coefficient (e.g., Ghisellini & Svensson 1991). 

In the SSC scenario, the relativistic electrons interact with synchrotron radiation pho- 
tons through the IC scattering. The IC emissivity is calculated by 

max 2 7 

^c(-c) = f /" ^/ -|^iV(7)/K.c)^W, (5) 

where Ui is the frequency of the incident photons emitted by the synchrotron radiation 
between i/™™ and u™^, /3 — v/c, and 72 are the lower and upper hmits of the scattering 
electrons, and /(i^i, I'c) is the spectrum produced by scattering monochromatic photons of 
frequency z/j with a single electron. The medium is transparent for the IC radiation field, 
hence /c(t'c) = e^i^^jR- 

Assuming that /g.c is an isotropic radiation field, the monochromatic luminosity around 
the source is obtained by 

L{5u) = 47r'R'hM5\ (6) 
Then the observed fiux density is given by 

47r^i?^/.,e(z/)^^(l + ^) 
= , (7) 

where o^l is the luminosity distance of the source and fobs — i^S/{l + z). 

In the GeV-TeV regime, the KN effect could be significant and make the IC spectrum 
have a high energy cut off. We take this effect into account in our model calculations. The 
high energy gamma-ray photons can also be absorbed by EBL, yielding electron-positron 
pairs, and the observed spectrum in the VHE band must be steeper than the intrinsic one. 
The absorption in the GcV-TeV band is considered with the EBL model as proposed by 
Franceschini et al. (2008). 



-6- 



3.2. Our Strategy and Results 

As shown above, the spectrum derived by the syn+SSC model is described with the 
following parameters, i.e., the strength of magnetic field B, the size of the radiating region 
R, the Doppler boosting factor 6, and the electron energy spectrum {pi, p2, 7min, 7b, 7max, 
and the electron density parameter Nq) (Tavecchio et al. 1998). Since S = [T{1 — /3cos9)]~^, 
where T = 1/^^1-/3'^ and I3^v/c, we take 5 ~ F in this paper because the relativistic jets 
of blazars are close to the line of sight (cos^ ~ 0). Besides pi and p2, these parameters are 
not directly observable. They are coupled with each other in the framework of the syn+SSC 
model. Our strategy and procedure to fit the SEDs are as follows: 

• Derive pi and p2 values from the spectral indices of the observed SEDs, which are 
obtained by fitting the observed data below and above the first bump with a power- 
law function. According to the syn+SSC model, the spectral indices for the two bumps 
are consistent. Therefore, if no data are available for the first bump, the observed data 
in the second bump is used. For six SEDs, the observed data are not enough to 
constrain pi, we search for historical pi from the published paper as listed in Table 1. 

• Set the values of 7inin and 7max- 7min is quite uncertain and is constrained with the 
observed SEDs via a method reported by Tavecchio et al. (2000). If no observational 
data are available to confine 7min, we take 7inin = 2 in our calculation. The 7inax can 
be constrained by the last observation point in the first bump of the SEDs, or be taken 
as 7inax = 1007b in this paper. 

• Set the radius (R) in the co-moving frame with the observed minimum variability 
timescale (At), i.e., R = c5At/{l + z). The values of At are listed in table 1 and 
the references are presented in Appendix. Rapid variabilities at timescales from hours 
to intra-day are observed in BL Lacs objects (Xie et al. 2001; Fossati et al. 2008), 
therefore, for some sources At is taken as one day when no variability timescales are 
available (see details in Appendix). 

• Set the lower limit of 5 with the 7-ray transparency condition for pair-production 
absorption (Dondi & Ghisellini 1995), 

S>Sr-[^A(^ + zr^]'/^''-''\ (8) 

where F[i>q) is the fiux density at the target photon frequency vq = 1.6 x 10^°/i/-y Hz, j3 
is the spectral index of the target photons (/3 = «i for uq < Ug and /3 = q;2 for uq > J^s), 
and tvar is the minimum variation timescale that is taken as tvar = At. 
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• Replace the values of 7b and No with the functions of B, S, the peak frequency z/g and 
the corresponding synchrotron radiation luminosity Lgy^ (e.g., see equations (2), (5) in 
Zhang et al. 2009) i.e., a set of B and 6 should reproduce the synchrotron radiation 
bump in a given parameters set of pi, p2, and R. 

• Constrain B and 6 with the ^^-minimization technique. Note that no error is available 
for some data points. We estimate the errors of these data points with the average 
relative errors of the data points whose errors are available. Since the relative errors 
of the data for the synchrotron bump are usually much smaller than that for the SSC 
bump, we derive the average relative errors of the data in the two bumps separately. 
We fit the distributions of the observed relative errors with a Gaussian function and 
take the central values of the Gaussian fits as the average relative errors. In case of a 
bow-tie in the SEDs, we use the average flux and its error at the low and high frequency 
ends as well as at the middle frequency of the bow-tie. We vary the values of 6 and B 
in wide ranges and calculate the corresponding of each parameter set of B and 6. 
Then, we derive the probability of the fit by p oc e~^'^^'^. We make the contours of p in 
the B — 6 plane and find the best fitting parameters. 

Following the above procedure, we finally get the contours of p for each SED in our sample. 
We find that the contours are usually thin and slanted ellipses, even degenerate into an 
unclosed ellipse for some SEDs that their SSC peaks are poorly observed. The contours of 
the SEDs in our sample are classified into three groups as following: 



• Closed ellipses in 1 a significance level. In this case, both B and S are constrained 
since the synchrotron radiation and SSC process peaks of the SED are well confined 
by the observation data. We get 16 such cases in our SED sample. As an example, 
the contours of p for the SED of 3C 66A are shown in Fig. [2](a). Both B and S can 
be constrained in 1 a significance level. We fit the the distributions of p for B and 5 
with a Gaussian function and derive the value of p that corresponds to 1 a significance 
level, as shown in Fig. [2](b), (c). We take the central values of the distributions of p 
for B and 6 as the centers of the contours and report our results in the form as df"^ 
and -B^"^ , where "up" and "low" mark the values at high and low ends of the contours 
in la significance level. The lower limit of 5 constrained with the 7-ray transparency 
condition for pair-production absorption is also displayed. Therefore, the lower limit 
of 5 should be taken as 5io„ = max[5T-, 5\o^] with the corresponding B^^. 

• Upward-opened ellipses in 1 cr significance level. In this case, the distributions of p for 
B and 5 still show as a hump that can be roughly fitted with a Gaussian function. The 
high end of 5 and the corresponding low end of B ((5up and -Blow) can be constrained 
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with the Icr contour. The low end of 5 and the corresponding high end of B cannot 
be constrained with the contour, but their hmits {6r and 5^) can be obtained with 
the 7-ray transparency condition. We get eight such cases in our SED sample. As an 
example, we show the contours of Mkn 180 in Fig. [2](d). We take the parameters of 
6 and B that present the minimum as the central values and report our results as 
6fr and 5<f- . 

• Downward-opened or unclosed ellipses. In this case, B and S cannot be confined by 
the contours. There are ten such SEDs in this group. As an example, we show the 
contours of Mkn 421 in the high state in Fig. EJ^g). Inspecting the distribution of p 
for S, one can find that p increases and comes to almost a constant at a given 6c- The 
constant value p roughly corresponds to the best fit to the data with the minimum 6. 
Therefore, we fit the distribution of p with a single-sided Gaussian function to obtain 
6c- Similarly, we get Be from the distribution of pioi B- We also derive the lower limit 
of 6 and the corresponding upper limit of B with the 7-ray transparency condition. 
Therefore, the limits of 6 and B are taken as 5 > max[5T-, 6c\ and B < min[i?^, B^]- 



As shown in Fig. [21 The contours of p are usually thin and slanted ellipses. The slopes 
of the major axis of the ellipses are usually -2 to -3, indicating that B and 6 are dependent 
to each other, i.e., B cx This relation seems to be consistent with that derived from 

the total luminosities of the two bumps (Lsyn and Lssc! Tavecchio et al. 1998)E, 



2L. 



1/2 



in the Thomson regime, and 

''g{ai,a2) 



^ Q ^ ^a^ 



(10) 



in the KN regime, where ^(ai, as) = exp[^ + ^^^^]- 

Our analysis shows that the major-axis of the contours is shorter when the sampling of 
a SED is better, hence tighter constraints on B and 6 are obtained, as shown in Fig. [2]^a) 
for 3C 66A. The best fitting parameters (or the limits) of B and 6 with the values of are 
listed in table 1 and the model curves derived with the parameters are also presented in Fig. 



^Note that the relation between B and S is not consistent with that derived from the two peak frequencies 
under the mono- frequency approximation condition, neither for Thomson regime nor for KN regime (e.g., 
see equations (4)(16) in Tavecchio et ah 1998). 
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[U We find eight SEDs in our sample that are also modeled with the single-zone syn+SSC 
model by other authors. We compare the derived B and 6 with that reported in literatures 
for these SEDs as following: 



• Mkn 421 in the high state 

BLazejowski et al. (2005) reported a parameter set of B = 0.26 G, 5 = 14, and R = 
7 X 10^^ cm (corresponding to At = 4.8 h), but they suggested that these parameters 
cannot explain the TeV data well. Therefore, they used a multi-zone SSC model to fit 
the SED. However, by our strategy we find that the SED can be well described using 
the single-zone SSC model with B = 0.14 G,S = 29,R = 9.1 x 10^^ cm, and At = 3 h. 
The most significant difference between our results with that reported by BLazejowski 
et al. (2005) is that the derived 6 with our strategy is almost twice of their value, 
whereas B is almost half of their value. 

• Mkn 501 in the high state 

Tavecchio et al. (2001) reported B = 0.32 G, 6 
sponding to At = 1.8 h). We get B = OA G, 6 -- 
h. Our results are roughly consistent with theirs 

• W Com in the high state 

Acciari et al. (2009) reported that B = 0.24 G, 6 = 20, and = 3 x 10^^ cm. We 
get B = 0.17 G, 6 = U, R = 1.6 x 10^^ cm, and At = 12 h. The difference between 
our results and theirs is the size of the radiating region. They set i? = 3 x 10^^ cm 
(corresponding to At = 1.5 h for S = 20). The strong variabilities on timescales of less 
than one day are observed at X-rays band, so we take At = 12 h in our calculation. 

• lES 1959-1-650 

It was reported that a parameter set of i? = 0.04 G, 6 = 20, and R = 5.8 x 10^^ 
cm (corresponding to At = 2.8 h) can well represent the SED of the object in the 
high state (Krawczynski et al. 2004) and a parameter set of 5 = 0.25 G, 6 = 18, 
R = 7.3 X 10^^ cm (corresponding to At = 4 h; Tagliaferri et al. 2008) for the low 
state. One can find that the reported magnetic field strength by these groups for the 
source in the low state is 6 times of that in the high state. We do not find such a 
significant difference for our results. We get B = 0.25 G, 5 = 12, i? = 1.2 x 10^^ cm, 
and At = 10 h for the high state and B = 1.1 G, 5 = 11, R = 1.1 x 10^^ cm, and 
At = 10 h for the low state. Please note that we take At = 10 h in our calculation, 
which is much larger that used by these authors. 



= 10, and i? = 1.9 X 10^^ cm (corre- 
= 15, R = 1.5 X 10^5 cm, and At = 1 
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• PKS 2155-304 in the low state 

Aharonian et al. (2009b) reported B = 0.018 G, 5 = 32, and = 1.5 x 10^^ cm 
(corresponding to At = 48.4 h), and our results are B = 0.16 G, 5 = 50, i? = 9.7 x 10^^ 
cm, and At = 2 h. The main difference between our results with theirs is that the 
adopted Ai in our calculation is much smaller than that derived with their fitting 
parameters. 

• 3C 66A 

Abdo et al. (2011) reported a much lower B and a corresponding larger 5 than the 
typical values of B and 5 for AGNs. They got B = 0.01 G, 5 = 50, = 1.1 x 10^^ 
cm (corresponding to Ai = 29.4 h). We get S = 0.2 G, 5 = 24 by taking Ai = 12 h 
(corresponding to i? = 2.2 x 10^^ cm in the co-moving frame). The derived R from our 
calculation is only 1/5 of that used by Abdo et al. (2011). 

• PKS 1424+240 

Acciari et al. (2010c) reported that 5 = 0.18 G, 5 = 45, and = 4.5 x 10^^ cm 
(corresponding to At = 13.9 h). We get B = 0.23 G, 5 = 33, i? = 5.7 x 10^^ cm, 
and Af = 24 h. Our results are roughly consistent with that reported by Acciari et al. 
(2010c). 

According to the comparison discussed above and our analysis results, one can find that 
a given observed SED usually can be represented by the different parameter sets. They are 

roughly consistent in sense of a larger 5 corresponding to a smaller B. However, the values of 
B and 5 from literatures do not exactly follow the relation oi B (x ^ since the adopted R 

(and the corresponding At) is different by different groups. The size of a radiating region is 
important in SED modeling. The previous authors usually took a given value of R to fit the 
observed SEDs. In our analysis we incorporate the observed minimum variability timescale 
and the value of 5 to constrain the size of radiation regions. Most importantly, we determine 
the parameter space for a SED by the minimization technique. To estimate the errors 
of the bolometric luminosity (I/boi), jet power (Pjet), and other physical parameters for our 
analysis below, we also derive the ranges of these parameters with the parameter sets of B 
and 5 that can represent the SED fits within 1 a range, and then derive the errors of these 
parameters and present in the following figures. Note that we do not consider the error bars 
in our correlation analysis since the data points usually have a large error and even only a 
limit is available. 
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4. Physical Properties of the Radiation Region 

4.1. Distributions of Model Parameters 

We show the distributions of B, S, 7;,, and R derived from the SEDs in our sample in 
Fig. [3], except for those that only limits are obtained from the observed SEDs. It is found 
that the values of B are clustered at 0.1 ~ 0.6 G and the range of 6 is from 5 to 50. The 
values of 7b vary from 10^ to 10^. The radius R for most sources are (4 ~ 40) x 10^^ cm. 
We also show the distributions of and in Fig. ^e). It is found that z/g ranges from the 
infrared to X-ray band (from 10^'^ Hz to 10^° Hz). The value of u^. covers a range from 10^° 
Hz to 10^^ Hz, but most of them are narrowly clustered at 10^^ ~ 10^^ Hz. 

The magnetic field of the radiation region is critical to understand radiation mechanisms 
of the sources. However, the origin of the magnetic field for the core regions is still uncertain. 
The narrow distribution of B might indicate that the magnetic field may be independent of 
the other jet properties. It is generally believed that the electrons are accelerated by the 
violent shocks via the Fermi acceleration mechanism, and then the value of 7b refiects the 
properties of the shock to some extent. No correlation between 7b and B is found. Therefore, 
the magnetic field of the core region may not be produced through the amplification of the 
interstellar medium magnetic field by the shocks. It is possible that the magnetic field comes 
from the vicinity around the central BH or from the accretion disk (Harris & Krawczynski 
2006). The strength of the magnetic field in the inner accretion disk region can be estimated 
by comparing it with that of BH binaries (e.g., Zhang et al. 2000), e.g., -Bagn /-^Binary oc 
(^AGN /^Binary )^ (10 cV/l kcV)^, whcrc T is the typical peak temperature of the accretion 
disk, and then -Bagn is ~ lO^'^i^Binary ~ 10"^ G. This is much larger than what we found in 
the jets. It is thus plausible that the magnetic fields in the jets are carried from the accretion 
flow, but diluted signiflcantly as the jets propagate and expand outwards. This may also 
explain the fact that the magnetic fields further out, e.g., in the jet knots and hot spots, are 
even much weakeJf]. 



4.2. Jet Power 

Jet power (Pjct) is essential to understand the physics of jet production. It can be 
estimated with Pjet = ttR^T'^c{U'^ + f/p + f/^), where F is the bulk Lorentz factor of the 
radiating region and U'^, {i = e, p, B) are the energy densities associated with the emitting 



•^Zhang et al. (2010) proposed that the magnetic field in the jet knots and hot spots may be due to 
amplification of the interstellar medium magnetic field by the violent shocks. 
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electrons U'^, the cold protons U'^ and magnetic field U'^ measured in the comoving frame 
(Ghisellini et al. 2010), which are given by 



U, = m,c' j N{^)^d^, (11) 

Ul = m^c' j N{^)dj, (12) 

U'^ = B'/St,, (13) 

assuming that there is one proton per radiating electron. These powers can be calculated with 
the SED fitting parameters. The radiation power is estimated with the observed luminosity, 



= nR^T^cU: = l' ^ ^, (14) 



4^4 - 4^2 ' 



where we take Lobs ~ -^boi (bolometric luminosity) and F ~ 5. may serve as a robust 
lower limit of Pjet (Celotti & Ghisellini 2008). We calculate Lboi in the band 10" - lO^^^ Hz 
based on our SED fits, then derive Pr- The calculated Lboi and Pjet are reported in Table 1 
and no statistical correlation between them is found (as shown in Fig. ID^a)). As the above 
statement, Pjet oc 5^, so the lower limit of 6 result in the lower limit of Pjet- The limits of 
Pjet for those sources, the SEDs of which can only provide a limit for 6, are also marked 
in the correlation figures, but are not considered in the following distribution figures. The 
distributions of the powers and their ratios to the total power, i.e., ee = Pe/Pjet, Cp = Pp/Pjet, 
= Pe/Pjet, and er = Pr/Pjet, are shown in Fig. |5l One can observe that the power in 
the jets is carried by cold protons for most sources, with a ratio = Pp/Pjet > 0.5 for most 
sources. The portions of the power carried by the magnetic field are comparable with e^.. 

It is generally believed that the jets are fed by accretion of the central massive BHs. 
With our results, we investigate the relation of Pjet to the mass of the central BH with 
a sub-sample of 18 sources in our sample. The masses of the BHs of these sources are 
collected from literature and reported in table 1. The jet power Pjet as a function of BH 
mass Mbh is shown in Fig. lU^b). A tentative anti-correlation between the two parameters is 
observed, i.e., Pjet oc Mj^^, a = 1.22 for high state and a = 1.21 for low state, respectively. 
The correlation coefficients estimated with the Spearman correlation analysis method are 
r = —0.74 and p = 0.04 for high state data and r = —0.86 and p = 0.007 for low state data. 
If the jet is powered purely by accretion of BH and refiects the accretion rate to some extent, 
i.e., Pjet = TjMc^, the negative correlation of jet power and BH mass would imply that the 
accretion rate is correlated with jj^- That is clearly unreasonable, and thus the relationship 
between Pjet and Mbh disfavors the pure accretion powered jet scenario. Therefore, we 
suggest that drawing the spin energy of the central BH should play a significant role in 
producing these jets. On the other hand. Fig. Ht^b) would imply a decrease in the BH spin 
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with an increase in Mbh- These suggest that the growth of BH mass in our sample (all with 
small redshift z < 1 ) is mostly through a series of random accretion processes that decrease 
the BH spin statistically (King & Pringle 2006; King et al. 2008, 2011). Wang et al. (2009) 
and Li et al. (2010) also suggested that minor mergers are important in the BH growth at 
low redshift and major mergers may dominate at high redshift, consistent with our results. 

However, for either the accretion or the spin model of BH to drive the jets, both of 
them face the cr-problem as reported by Harris & Krawczynski (2006), where a is the ratio 
of electromagnetic energy density to particle energy density. Both models require that the 
jet formation needs a Poynting flux dominated energy transport by a strongly magnetized 
or high-(j plasma; however, parsec-scale observations indicate that the jets are dominated 
by particles and are low-a plasma. We also present the correlation between and Pb 
(excluding the data with upper or lower limits) as shown in Fig. 111(c), and find that the jets 
of those sources are indeed dominated by the particles too, either for the electron pair model 
or for the electron-proton model. Our results thus also face the cr-problem. As described 
in section 4.1, the magnetic field in the jets may be carried from the accretion flow (Harris 
& Krawczynski 2006) and are diluted significantly as the jets propagate. It is therefore 
required that the majority of the Poynting flux is converted to particle energies efficiently 
even beyond equipartition. However it remains unclear how and why this conversion is so 
efficient. 



4.3. Flux Variation 

The SEDs of ten objects in our sample are obtained in the low and high TeV states. We 
compare Us and u^. between the high and low states in Fig. El It is clear that the SEDs in the 
high states shift to a higher energy band. Therefore, the flux variation is usually accompanied 
by a spectral shift. To investigate what may be responsible for the flux variations and the 
spectral shift in the low and high states, we derive the ratios of the flux density at 1 TeV 
{Ri Tev) and the physical parameters [Rx) in the high state to that in the low state for the 
ten sources, where x stands for Lboi, B, 6, 7b, and Pjet- Fig. [71 shows Rx as a function of 
Ri TeV- Note that 7^ = 3 7x^06 'm ^^'^ ^ ^ so the lower limit of 6 and the upper 

limit of B result in a lower limit of 7^ and the lower/upper limits of the ratios of 7?, for 
high to low states. The limits of Rx resulting from the limits of S and B are also marked 
in the following correlation flgures. Here, the limits of Rx are also derived by the contours 
of B and 5; taking the lower limit of Rs for Mkn 421 as an example, we get Rs = 29/43. 
The values of 7b, -Z^bob and Pjet of the high states are systematically higher than that of 
the low states. We measure the correlations between Ri xev and Rx with the Spearman 



- 14 - 



correlation analysis method, and find a tentative correlation between Ri xev and R^^^, with 
the correlation coefficient r = 0.66 and the chance probability p = 0.04. The best fit line 
between R^^ and i?i xev is also shown in Fig. EJ^e). No statistically significant correlations 
between Ri xev and other parameters are found. It is well known that the electrons are 
accelerated by the violent shocks via the Fermi acceleration mechanism, and then the value 
of 7b may reflect the different microphysical properties of the shock to some extent. As 
shown in Fig. [3]^c), the range of 7b is from 10^ to 10^ and the values of 7b among sources, 
even among states of a given source, are different, likely indicating that the properties of 
shocks significantly vary among sources (or states). Therefore, it is possible that the flux 
variation is probably due to the variation of the physical properties for the shocks, such as 
the obliquity of the shock and the level of magnetic field turbulence in the shock (Summerlin 
& Baring 2012). An interesting phenomenon for the flux variation of an individual source 
is that the increase of 7b is accompanied with the decrease of B as shown in Fig. Ulh). A 
detailed analysis on this phenomenon will be presented in our next paper. 



5. Summary 

The peak frequency and the peak luminosity of the high energy bump of SEDs for 
AGNs can be well determined with simultaneous or quasi-simultaneous observations in the 
TeV and GeV bands. We compile the broadband SEDs of 24 BL Lac objects that were 
simultaneously or quasi- simultaneously observed with Fermi/LAT and other instruments 
from literature. The clean SEDs without contaminations of the accretion disks and external 
IC processes of these sources are good candidates for investigating the radiation mechanisms 
and the physical properties of the AGN jets. Our results are summarized as follows: 

• We find that the one-zone syn-|-SSC model can well represent the observed SEDs. No 
excess due to the interaction between the TeV photons and the extragalactic back- 
ground light in GeV band over SSC component is found for the sources in our sample. 

• The distribution of 7b ranges from 10^ to 10^, but the magnetic field strength B is 
clustered within a narrow range of 0.1-0.6 G. These results indicate that the properties 
of shocks significantly vary among sources (or states), and the magnetic field may not 
be due to the amplification of the interstellar magnetic field by the shocks in the jet. 
We propose that the magnetic field of the core region may originate from the accretion 
flow. 

• The Doppler boosting factor 6 of the jets ranges from 5 to 50, and the sizes of the 
radiating regions are roughly (4 ~ 40) x 10^^ cm. Significant flux variations are observed 
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for the sources in our sample. The SEDs in the high state shift to higher frequencies. 
A tentative correlation between the ratio of the flux density at 1 TeV and the ratio of 
the 7b in the low and high states is found, indicating that the flux variation is probably 
due to the variation of the physical properties for the shocks, such as the obliquity of 
the shock and the level of magnetic field turbulence in the shock. 

• We calculated the bolometric luminosity and the jet power for the sources in our 
sample. The jet power is dominated by the kinetic energy for most sources. No 
correlation between Lboi and Pjet is found. An anti-correlation between the jet power 
and the mass of the central BH is observed, i.e., Pjet oc Mg^. This disfavors the scenario 
of a pure accretion-driven jet. We suggest that extraction the rotation energy of the 
central BHs would be significant for these sources, and that thus BHs with smaller 
masses may be spinning more rapidly in this sample. 

6. Discussion 

With a large sample of different types of blazars, Fossati et al. (1998) reported a 
spectral sequence of FSRQ-LBL-HBL, i.e., along with this sequence, an increase in the 
peak frequencies corresponds to the decreases of bolometric luminosity and the ratio of the 
luminosities for the high- and low-energy spectral components. This spectral sequence was 
interpreted by Ghisellini et al. (1998) with the evolution of the external photon fields (such 
as the BLRs) of the blazars, or more physically, the sequence is due to the different mass 
and accretion rate of the BH (Ghisellini & Tavecchio (2008). More recently, Chen & Bai 
(2011) extended this sequence to narrow line Seyfert 1 galaxies. The significant cooling 
effect for the electrons by the external field photons from the BLR of FSRQs may result in 
a low 7b in the electron spectrum. The contribution to the luminosity from the IC process 
of the external photon fields would significantly increase the total luminosities and the ratio 
of Lic/Lsyn- Therefore, observationally, one can expect that both Lboi and Ljc/i^syn are 
anti-correlated with i/g, and thus an anti-correlation between 7b and U[ in the co- moving 
frame is also expected (Ghisellini et al. 1998). Lboi and Lic/Lsyn as a function of Us are 
shown in Fig. |8l 7b as a function of Ul {U[ = f/syn,avaii + in this work, where f/gyn^avaii 
the available photon energy density of synchrotron radiation for IC process) is also shown in 
Fig. [HI No correlation is found for pairs Lboi — i^s and 7b — U[. This lack of correlation in our 
sample can be explained by the fact that the EC process is negligible for the BL Lac objects 
in our sample. We also show the sample of Ghisellini et al. (2010) in the 7b-t^r plane (gray 
circles in Fig. El^a)). Although the correlation between 7b and U[ of our sample sources has 
large scatters, comparing our results with that of Ghisellini et al. (2010), the two results are 
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consistent and the sources in our sample distribute in the left top of Fig. [Uj^a), where the 
EC process is not important. 

However, the ratio of Ljc/i^syn is indeed anti-correlated with z/g in their high and low 
states, except for PKS 2005-489 in the high state. The Spearman correlation analysis yields 
a correlation coefficient r = —0.59 and a chance probability p = 0.07 for the high state data, 
r = —0.9 and p = 3.4 x 10~^ for low state data, respectively. Excluding PKS 2005-489 
in the high state, our best fits give log Lic/i^syn = (3.27 ± 1.) — (0.23 ± 0.06)logz/s and 
log Lic/i^syn = (3.5 ± 0.95) — (0.2 ± 0.06) logi/g for the sources in the low and high states, 
respectively. Because all the sources in our sample are BL Lac objects, the external photon 
fields outside their jets are much weaker than the synchrotron radiation photon field and the 
EC process is thus not considered in this work. Therefore, the Lic/-Z^syn ~ anti-correlation 
may have a different physical origin from the blazar sequence. It is possible that this anti- 
correlation is due to the KN effect. As increases, the SEDs shift to the higher frequency 
end and the KN effect should be more significant. Thus, the ratio of Lic/i^syn would decrease 
as Us increases since Ub is almost constant. 

The TeV photons from the high redshift universe can be absorbed through interac- 
tion with the extragalactic background light (EBL), producing electron-positron pairs. The 
electron-positron pairs may scatter the CMB photons to the GeV band. Therefore, the in- 
trinsic spectrum at the TeV band should be harder than the observed one, especially for the 
sources at high redshift, and an excess component in the GeV band may be observed in the 
observed SEDs as suggested by some authorfl The power of such a component depends on 
the EBL modei^l and the intergalactic magnetic field (IGMF). Assuming that the intrinsic 
spectrum of TeV emission is the same as that of the GeV emission, one can confine the EBL 
model and the IGMF. Using the data of Mkn 501 during its high state, Dai et al. (2002) 
calculated the spectrum under different IGMF strengths and reported that the cascade GeV 
emission is detectable by Fermi/LAT if the IGMF strength is < 10^^^ G. Comparing the 
predicted emission in the GeV band with the upper limits measured by the Fermi/LAT for 
the blazar lES 0229-^200, Tavecchio et al. (2010a) suggested that the IGMF is larger than 
10~^^G. Our systematical analysis of the TeV BL Lac objects presented here clearly indicates 
that no excess emission over the SSC component is observed: this may be due to a stronger 



''it was also reported that the electron-positron pairs would be deflected by the intergalactic magnetic 
field (IGMF) from the initial TeV photon direction and thus the secondary GeV photons produced by them 
would exhibit a halo around the central bright point-source instead of a GeV excess in the SEDs (Aharonian 
et al. 1994; Dolag et al. 2009; Neronov & Semikoz 2009). 

■^Using the multi-wavelength galaxy observation data and the luminosity functions in the near-IR band, 
Dominguez et al. (2011) derived the EBL model to redshift 4. 
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IGMF than previous predictions (Dai et al. 2002; Yang et al. 2008; Tavecchio et al. 2010a). 

7. Appendix 

Mkn 421. It is the first TeV AGN confirmed by Whipple (Punch et al. 1992) and is 
a high-frequency-peaked BL Lac (HBL) object. Violent variation of the fiux at GeV-TeV 
regime was detected and associated with that observed at the X-ray band (BLazejowski et 
al. 2005). The data of high state are from BLazejowski et al. (2005), when the source 
underwent an outburst in 2004 April with the peak fiux ~ 135 mCrab in the X-ray band and 
~ 3 Crab in the gamma ray band. The variabilities on timescales of hours were observed, so 
we take At = 3 h to model the SED of the high state (to see Figure 10 in BLazejowski et 
al. 2005). The data of the low state are taken from Paneque et al. (2009), when Mkn 421 
is found in a rather low/quiet state from August 2008 to August 2009. At for modeling the 
low state SED is also taken as 3 h. The high energy gamma-ray emission of this source has 
also been detected by the ARGO-YBJ detector (Aielh et al. 2010, Bartoli et al. 2011). 

Mkn 501. It was detected by Whipple during an observation of 66 h. An average 
fiux of (8.1 ± 1.4) X 10~^^cm~^s~^ above 300 GeV and the variability on timescales of days 
were observed (Quinn et al. 1996). From 1997 April to 1999 June, the observations with 
BeppoSAX showed that the peak frequency of the synchrotron emission shifted from 100 
keV back to 0.5 keV, and correspondingly the fiux decreased (Tavecchio et al. 2001). A 
multi-wavelength campaign in 2006 July with Suzaku and MAGIC was performed. During 
this program the average VHE fiux above 200 GeV is ~ 20% of Crab fiux with a photon index 
2.8 ±0.1 from 80 GeV to 2 TeV, indicating that the source was in a low state (Anderhub 
et al. 2009a). At is taken as 3 h (Figure 2 in Anderhub et al. 2009a) to fit this observed 
broadband SED in this paper. The data of the high state SED (with At=l h) are from 
Tavecchio et al. (2001). The data of Fermz/LAT observations are from Abdo et al. (2009), 
but they are not simultaneous with the data of low and high states. 

W Com. The first TeV intermediate-frequency-peaked BL Lac (IBL) object, an evi- 
dently bright outburst in the optical and X-ray bands was observed in 1998 (Tagliaferri et 
al. 2000). W Com is confirmed to be a TeV source by the observations of VERITAS during 
a strong TeV fiare in the middle of March 2008 with an integrated photon fiux above 200 
GeV of ~ 9% Crab (Acciari et al. 2008) and the data of the TeV observation are presented 
as opened squares in Fig. [1] Both the SEDs quasi-simultaneously obtained during the TeV 
fiare and during the optical/X-ray outburst can be well fitted by the single- zone syn+SSC 
model (Zhang 2009). The broadband SED quasi-simultaneously obtained with observations 
of Fermi/L AT and swift from Tavecchio et al. (2010b) is also shown in Fig. [T](red squares). 
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Subsequently, another outburst of very high energy gamma-ray emission was detected in 
2008 June by VERITAS with the flux of (5.7 ± 0.6) x lO-^^cm-^s-^ (blue squares in Fig. 
[H Acciari et al. 2009a), three times brighter than the observation of March 2008. The flux 
upper-limit derived by quasi- simultaneous observation of AGILE is also presented in Fig. [T] 
(gray star). The strong variabilities on timescales of less than one day are observed in X-rays 
band (Fig. [T]in Acciari et al. 2009a), so we take At = 12 h in this paper. For the SEDs of 
high state, Acciari et al. (2009a) pointed out that both the SSC model and the SSC-I-EC 
model can explain the observed SED well, but the SSC+EC model makes the magnetic field 
strength more close to the equipartition condition. 

BL Lacertae. The first confirmed TeV low-frequency-peaked BL Lac (LBL) object, its 
VHE gamma-ray emission was detected by MAGIC during 2005 August to December with 
an integral flux of (0.6 ± 0.2) x lO^^^cm^^s^^, corresponding to 3% Crab flux (Albert et al. 
2007a). The photon index from 150 to 900 GeV is rather steep with F = —3.6 ± 0.5, and 
the light curve shows no significant variability. The simultaneous observation for MAGIC 
in optical band was performed by KVA. The broadband SED together with EGRET data 
(opened squares) in 1995 is shown in Fig. [T] (blue squares). During an optical outburst in 
July 1997, BL Lacertae was detected in X-ray band by RXTE and in gamma-ray band by 
EGRET, implying that the source was bright and variable at the both bands (Madejski et 
al. 1999). The spectra in X-ray and gamma-ray bands are hard. Madejski et al. (1999) 
reported that the X-rays are produced by SSC process while the gamma-rays are produced 
by Comptonization of the broad emission line flux. They also suggested that the the emitting 
regions correspond to the collision of two blobs or shells moving at slightly different velocities, 
therefore, the different SEDs correspond to different collisions, which can occur at different 
location in the jet. If the emission regions are inside of the BLRs, they would emit most of 
their energy in the GeV band. The multiwavelength data during Fermi/LAT observation 
(red squares in Fig. [H Tavecchio et al. 2010b) are also considered, but no simultaneous 
data in TeV regime are obtained. A rapid variability was detected in the soft X-ray band 
on timescales of 3-4 hours (Ravasio et al. 2002), therefore, we take At = 3 h to fit the 
broadband SEDs of the two states. 

PKS 2005-489. A HBL object was detected by HESS in 2003 and 2004 (Aharonian et 
al. 2005). A significant signal of VHE emission was detected in 2004 with the integral flux 
above 200 GeV of ~ 6.9 x 10^^^ cm^^ s^^, corresponding to 2.5% Crab flux. The flux level 
in 2003 was lower than that measured in 2004, indicating that the activity in 2004 increased, 
but no significant variability on timescale less than a year was found. The multiwavelength 
observation by HESS, XMM-Newton and RXTE satellites from 2004 to 2007 was performed 
(HESS Collaboration 2010a), indicating that the large flux variations at the X-ray band are 
coupled with weak or no variations in the VHE band. A broadband SED with the low X-ray 
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flux is taken as the low state and shown as red symbols in Fig. [H A very high state of X-ray 
emission for this source in 1998 was detected by BeppoSAX and RossiXTE (Tagliaferri et al. 
2001). HESS and Fermi/LAT simultaneously observed the source in 2009, when the X-ray 
flux was comparable to the flux level in 1998 (Kaufmann et al. 2010). Therefore, we compile 
the data of the two observations together and take it as the high state (blue symbols in Fig. 
[1]). The variability timescale for modeling the two SEDs is taken as At = 10 h (Tagliaferri 
et al. 2001). 

lES 1959+650. A HBL object with a strong TeV outburst in May 2002 was detected 
by VERITAS (Holder et al. 2003), HEGRA (Horns et al. 2002) and CAT (Aharonian et 
al. 2003). During the outburst, the flux level reached to 3 times the Crab flux and there 
was evidence for strong variability. The multiwavelength campaign in radio, optical. X-ray 
and TeV gamma-ray bands was performed from 2002 May 18 to August 14 (Krawczynski 
et al. 2004). A time-averaged spectrum corresponding to the TeV gamma-ray high state is 
considered in this work (blue squares in Fig. [1]). The observations show that X-ray flux and 
gamma-ray flux had tentative correlation, but no correlations of optical variability with X- 
ray and gamma-ray emission were found. Another multiwavelength campaign was performed 
in 2006 May, when the source exhibited a high state at optical and X-ray bands and was in 
the lowest level at the VHE band (Tagliaferri et al. 2008). There were variabilities in optical 
and X-ray bands. The spectrum butterfly of Fermi/LAT observation is also shown in Fig. 
[H The fastest gamma-ray flux variability timescale of 10 hours based on the Whipple data 
is obtained and then is taken as the value of At to model the SEDs of the two states. 

lES 2344+514- The third BL Lac object was detected at the VHE band with Whipple. 
The detection of VHE emission mostly came from an apparent flare on 1995 December 
and the average flux above 350 GeV was (6.6 ± 1.9) x 10^^^ cm~^ s^^, 63% of the Crab 
flux (Catanese et al. 1998). The observation of MAGIC between 2005 August 3 and 2006 
January 1 presented a steep spectrum with photon index F = —2.95 and a flux 6 times 
below the 1995 flare, indicating that the source was in low state (Albert et al. 2007b). 
No evidence for variability was found during the MAGIC observations. The simultaneous 
optical observation with MAGIC was performed by KVA. The broadband SED is shown as 
red solid squares in Fig. [TJ The data of Swift satellite observation (red opened squares) on 
19 April 2005 (Tramacere et al. 2007) and the bow-tie of Fermi/LAT observation (Abdo 
et al. 2009) are also presented to constrain the radiation model. A broadband SED for a 
high state quasi- simultaneously obtained by VERITAS and Swift on December 2007 is also 
considered in this work (blue squares), and the data are from Acciari et al. (2011). The 
measured flux above 300 GeV was (6.76 ± 0.62) x 10~^^ cm~^ s~^, corresponding to 48% 
of Crab flux. The highest X-ray emission ever observed was measured by Swift/XRT, and 
was correlated with the VHE gamma-ray emission. The TeV data of the high state in Fig. 
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[T]have been corrected for absorption by EBL. The variabilities on timescale of a few hours 
(Albert et al. 2007b) and on timescale of ~ day (Acciari et al. 2011) are observed, so we 
take At as 12 hours for modeling the SED of two states. 

PKS 2155-304. A HBL object, its strong VHE emission was first detected in 1997 
November, at the same time the strongest X-ray emission ever observed and GeV gamma 
rays were also detected by BeppoSAX and EGRET (Chadwick et al. 1999). A simultaneous 
observation with HESS, Chandra and the Bronberg optical observatory was carried in the 
night of July 26-30 2006 when the source was in a high-activity state at gamma-ray band 
and the gamma ray flux reached ~ 11 times the Crab flux (Aharonian et al. 2009a). The 
emission between X-ray and VHE gamma-ray is strongly correlated. The broadband SED of 
this campaign is shown as blue squares in Fig. [TJ The variability on timescale of 2 hours is 
observed at optical (and higher) frequencies (Aharonian et al. 2009a), thus we take At = 2 
h. The GeV-TeV observation with Fermi and HESS was performed between 25 August and 
6 September 2008, and the low-energy component was simultaneously covered by ATOM 
telescope, RXTE and Swift observations (red squares in Fig. [1], Aharonian et al. 2009b). 
During that period PKS 2155-304 was in a low-activity state at X-ray and gamma-ray bands, 
whereas the optical flux was much higher. The optical emission was correlated with the VHE 
emission, but no correlation between X-ray and VHE was found. At is also taken as 2 hours 
to model the low state SED. 

lES 1101-232. A HBL object hosted by an elliptical galaxy, it was detected by HESS 
in March- June 2004 and 2005 with a very hard spectrum and no significant variation was 
found (Aharonian et al. 2007a). A multiwavelength campaign data, including observations 
of VHE by HESS and X-ray by RXTE satellite in 2005 and XMM-Newton in 2004, were 
obtained (red squares in Fig. [H Costamante 2007). The Suzaku observation simultaneously 
covered with the HESS measurement was carried out in 2006 May (blue squares in Fig. [H 
Reimer et al. 2008). No significant X-ray or gamma-ray variability was detected during this 
program, and the object was in a quiescent state with the lowest X-ray flux ever measured. 
lES 1101-232 was not detected by Fermi/LAT (Tavecchio et al. 2010b), and only an upper- 
limit was given as shown in Fig. [H No variability timescale is found for the corresponding 
SEDs, so we take At = 12 h to fit the SEDs of two states. 

S5 0716-h714. It is a LBL object. The MAGIC observations were performed in Novem- 
ber 2007 and in April 2008 and detected its TeV emission (Teshima et al. 2008; Anderhub 
et al. 2009b). The integral flux above 400 GeV was ~ 7.5 x 10~^^ cm~^ s~^, corresponding 
to 9% Crab flux. The optical emission of S5 0716-1-714 was simultaneously observed by KVA 
when the source was in a high state at optical band, and most of the gamma-ray emission 
signal came from the phase of the optical high state of the object, suggesting a possible 
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correlation between the VHE emission and optical emission (Anderhub et al. 2009b). The 
data quasi- simultaneously obtained in April 2008 by KVA, Swift and MAGIC are considered 
and shown as blue symbols in Fig. [TJ The data of Fermi/'LP^I observation with quasi- 
simultaneous Swift observation from Tavecchio et al. (2010b) are also considered and shown 
as red symbols in this work. At for the high and low states is taken as 3 hours (Foschini et 
al. 2006). 

3C 66A. A IBL object, the observation of VHE was performed from September 2007 
through January 2008 by VERITAS and it was confirmed to be a TeV source with an integral 
flux above 200 GeV 6% of Crab flux (Acciari et al. 2009b). The observed spectrum can 
be characterized by a soft power-law with photon index F = 4.1 and a variability on the 
timescale of days was found. The simultaneous GeV- TeV observations by Fermi/LAT and 
VERITAS were performed in October 2008 (Reyes et al. 2009). The simultaneously observed 
broadband SED is from Abdo et al. (2011). The variability on timescales of less than a day 
is observed (Abdo et al. 2011) and then At is taken as 12 hours. Abdo et al. (2011) 
reported a SSC+EC model to describe the SED since the magnetic fields are significantly 
closer to equipartition than in the pure SSC case, i.e., Ls/Le ~ 0.11 for SSC+EC model 
and Lb/L^ ~ 1.1 x 10~^ for pure SSC model. 

PG 1553+113. Its VHE emission with a very soft spectrum (photon index of F = 
4.0 ±0.6) was detected in 2005 by HESS, and no evidence of variability was found (Aharonian 
et al. 2006a). The VHE emission was subsequently confirmed by MAGIC (Albert et al. 
2007c), and the integral flux levels were consistent with that of HESS. PG 1553+113 is in 
the Fermi LAT bright AGN source list, but was not detected by EGRET because it was in a 
low state during the observation. This source is a HBL object and a bright X-ray source with 
many observations, but no strong or fast variability was detected in the X-ray band (Reimer 
et al. 2008). The redshift of PG1553+113 is unknown and the VHE observation indicates 
that the redshift is greater than 0.25. In this work, we take z = 0.3. The observations by 
Fermi/LAT from 4 August 2008 to 22 February 2009 show that it was a steady source with a 
hard spectrum in Fermi/LAT energy band (Abdo et al. 2010c). The data of broadband SED 
(black squares) are from Abdo et al. (2010c), including quasi-simultaneous observations by 
KVA, Suzaku, MAGIC and HESS in July 2006 and the observation of Fermi/LAT. At is 
taken as one day. 

lES 1218+30. 4. A HBL object, it was confirmed to be a TeV source by MAGIC 
in 2005 January, but no variability on timescales of days was found within the statistical 
errors (Albert et al. 2006b). In optical band, KVA observed the source simultaneously with 
MAGIC. From 2008 December to 2009 May, VERITAS monitored the source and revealed 
a prominent flare reaching ~ 20% of the Crab flux (Acciari et al. 2010a). The light curve 
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of TeV emission for this source showed day-scale variabihty. Here, we take At — 12 h. The 
observational flux of VERITAS is comparable to MAGIC (Weidinger & Spanier 2010). Swift 
observed this source between March and December 2005, and was quasi-simultaneous with 
the observation of MAGIC (Riiger et al. 2010). In this work, the data of the broadband 
SED include the observations of VERITAS, Fermi/LAT (from Abdo et al. 2009), swift and 
KVA. 

lES 1011+496. A HBL object, it was observed by MAGIC from 2007 March to May 
after an optical outburst in March 2007, and an integral flux above 200 GeV of 1.58 ±0.32 x 
10~^^ cm~^ s^^ was obtained (Albert et al. 2007d). The broadband SED is composed of the 
quasi-simultaneous observations of Swift and Fermi/LAT (from Tavecchio et al. 2010b) and 
the observation of MAGIC in 2007 (opened squares) . No corresponding variability timescale 
for the SED is found, so we take At — 24 h. 

PKS 1424+240. A HBL object with an unknown redshift, its VHE emission was de- 
tected by VERITAS with a flux normahzation at 200 GeV of ~ 5.1 x 10"^^ TeV"^ cm-^ 
s~^ (Acciari et al. 2010c). The photon spectrum above 140 GeV can be described well by 
a power law with F ~ 3.8. During the period from February 2009 to June 2009, the flux 
of VHE emission was steady and the contemporaneous observation of Fermi/LAT also did 
not detect any variability (Acciari et al. 2010c). The broadband SED is established by si- 
multaneous observations of VERITAS, Fermi/LAT, Swift and MDM. Considering the EEL 
absorption, a redshift upper limit of 0.66 is inferred (Acciari et al. 2010c) and z = 0.5 is 
taken in this work. The X-ray variability timescale of about a day (Acciari et al. 2010c) is 
taken as the value of At to model the SED in this paper. 

RGB J0710+591. A well known HBL object, it was not detected by EGRET. It was 
observed at the VHE waveband by VERITAS between December 2008 and March 2009, and 
was conflrmed to be a TeV source (Ong et al. 2009, Acciari et al. 2010b). The observed 
spectrum from 0.31 to 4.6 TeV can be flt by a power law with a photon spectral index 
~ —2.69, and the integral flux above 300 GeV is 3.9 ± 0.8 x 10~^^ cm~^ s~^, corresponding 
to 3% of Grab flux. The VERITAS observation was complemented by contemporaneous 
observations from Fermi/LAT, Swift and Michigan-Dartmouth-MIT observatory (Acciari et 
al. 2010b). No corresponding variabihty timescale for the SED is found, so we take At — 24 
h. 

lES O8O6+524. A HBL object was detected by VERITAS in VHE gamma-ray regime 
between November 2006 and April 2008, and no signiflcant variabihty on months timescale 
was found (Acciari et al. 2009c). The observed photon spectrum from November 2007 to 
April 2008 can be fitted by a power law with F ~ —3.6 between 300 to 700 GeV. The integral 
flux above 300 GeV is ~ 2.2 x 10~^^ cm~^ s~^, corresponding to 1.8% of Crab flux. The data 
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obtained quasi-simultaneously by Swift and VERITAS observations (Acciari et al. 2009c) 
and the spectrum butterfly of Fermi/LAT observations (Abdo et al. 2009) are considered in 
this work. No corresponding variability timescale for the SED is found, so we take At — 24 
h. 

Mkn 180. A HBL object was detected with VHE gamma-ray emission by MAGIC 
during an optical outburst in 2006. The integral flux above 200 GeV is (2.3 ± 0.7) x 10~^^ 
cm~^ s~^ and corresponds to 11% Crab flux. The observed spectrum was rather soft with 
a photon index of F ~ 3.3 ± 0.7, and no variability was found (Albert et al. 2006a). A 
multi- wavelength campaign in 2008 was performed (Rugamer et al. 2011). The broadband 
SED is compiled with the observations of swift and Fermi (grey shaded) in 2008 and the 
observation of MAGIC in 2006. The bow-tie of Fermi/LAT observation for 1-year Catalog 
(dashed hue, Abdo et al. 2009) is also given in the SED. According to the light curve of 
Mkn 180 in Albert et al. (2006a), At is taken as one day in this paper. 

H 2356-309. A HBL object was detected with VHE emission by HESS from June to 
December 2004, with a integral flux above 200 GeV of 4.1 ±0.5 x 10~^^ cm~^ s"^ (Aharonian 
et al. 2006b). A simultaneous observation with HESS in 2004 at lower energy bands was 
performed by ROTSE-III (optical) and RXTE (X-rays). A broadband SED obtained simul- 
taneously by XMM-Newton and HESS on June 2005 (HESS Collaboration et al. 2010b) is 
presented in this work, in which the flux levels in X-ray and TeV bands are comparable with 
that measured in 2004. In the VHE band, signiflcant small-amplitude variations on time 
scales of months and years were detected. The X-ray measurements show that it was in a 
low state in this band. Unfortunately, the observations of Fermi/LAT only give an upper 
limit (Tavecchio et al. 2010b). At is taken as 24 h in this paper. 

lES 0347-121. A BL Lac object was detected by HESS between August and December 
2006 with an integral flux corresponding to 2% of Crab flux (Aharonian et al. 2007b). The 
photon spectrum from 250 GeV to 3 TeV can be described by a power law with a photon 
index F ~ 3.1. During the VHE observation, no significant variability was detected, so we 
take At — 12 h. The broadband SED was compiled with the data from the simultaneous 
observations of HESS, Swift and ATOM (Aharonian et al. 2007b). The upper limit detected 
by Fermi/LAT is from Tavecchio et al. (2010b). 

lES 0229+200 A HBL object was observed by HESS in 2005/2006, and confirmed to 
be a TeV source (Aharonian et al. 2007c). The integral fiux above 580 GeV is ~ 9.4 x 10"^^ 
cm~^ s~^, corresponding to ~ 1.8% of Crab fiux, and the observed spectrum is characterized 
by a hard power law with F ~ 2.5 from 500 GeV to 15 TeV. During the observation, no 
signiflcant variability on any scale was detected. Except for the data of the HESS observation, 
the broadband SED considered in this work also includes the data of Swift observation in 
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August 2008 and an upper limit of Fermi/LAT observations (Tavecchio et al. 2010b). No 
corresponding variability for the SED is found, so we take At — 24 h. 

RGB J0152+017. A HBL object was detected by HESS in late October and November 
2007 (Aharonian et al. 2008). The observed spectrum is well fit by a power law with 
r ~ 2.95, and the integral flux above 300 GeV corresponds to ~ 2% of Crab flux. The 

data of the broadband SED also include the simultaneous observations of Swift and RXTE 
(Aharonian et al. 2008), and the upper limit detected by Fermi/LAT (from Tavecchio et al. 
2010b). No corresponding variability for the SED is found, so we take At = 24 h. 

PKS 0548-322. A HBL object was observed between October 2004 and January 2008 
with the HESS, and confirmed to be a TeV source (Superina et al. 2008, Aharonian et 
al. 2010). The integral flux above 200 GeV is 1.3% of the Crab flux and the observed 
spectrum is characterized by a power-law with a photon index F ~ 2.86. Contemporaneous 
UV and X-ray observations in November 2006 were made by Swift, but it was not be detected 
by Fermi/LAT (Tavecchio et al. 2010b). No signiflcant variabihty was detected by HESS 
and Swift. In this work, the broadband SED includes the data of the Swift and HESS 
observations, together with the upper limit of Fermi/LAT observation. At is taken as 24 h. 

H 1426+428. A HBL object with the strongest TeV emission was detected in 2000 and 
2001 by Whipple with an integral flux of (2.04 ± 0.35) x 10"^^ cm-^ s"^ above 280 GeV 
(Horan et al. 2002). The object was monitored by Whipple from 1995 to 1998 during a 
general blazar survey, but no statistical signal was detected. No simultaneous broadband 
data are found and the data of broadband SED are from Wolter et al. (2008). The spectrum 
butterfly of the Fermi/LAT observations (Abdo et al. 2009) is also taken into account. No 
corresponding variability for the SED is found, so we take At — 24 h. 
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"^Thc state of each source in TeV band. "H" indicating "high state" and "L" indicating "low state" . ^ 



^z: rcdshift. 

''pi for six SEDs is cited from the references. Corresponding to: (1) Albert et al. (2007a); (2) Tavecchio et al. (2010b); (3) Woltcr ct al. (2008). 

''BH mass for 18 sources, among which fourteen are from Woo & Urry (2002), two are from Wu et al. (2002), and two arc from other references. 
The superscript "Re" denotes the reference given in column (16) and the superscript "Wu" denote the reference Wu et al. (2002). 

Note. — Columns: (4) (5) The energy indices of electrons below and above the break; (6) The minimum variability timescale; (7) The minimum 
Lorentz factor of electrons; (8) The beaming factors 5; (9) The inferred magnetic field strength; (10) oi SED fitting results; (11) The lower 
limit of 6 derived by the 7-ray transparency condition; (12) Flux density at 1 TcV; (13) The bolometric luminosity of each SED; (14) The total 
jet power; (15) The mass of the black hole in the center of each AGN; (16) The references. 

References. — (1) BLazejowski et al. 2005; Paneque et al. 2009; (2) Tavecchio et al. 2001; Anderhub et al. 2009a; (3) Tagliaferri et al. 2000; 
Acciari et al. 2008; Acciari et al. 2009a; Tavecchio et al. 2010b; (4) Albert et al. 2007a; Ravasio et al. 2002; Tavecchio et al. 2010b; (5) HESS 



Collaboration 2010a; Tagliaferri et al. 2001; Kaufmann et al. 2010; (6) Krawczynski et al. 2004; Tagliaferri et al. 2008; (7) Albert et al. 2007b; 
Tramacere et al. 2007; (8) Aharonian ct al. 2009a; Aharonian ct al. 2009b; Ricgcr & Volpc 2010; (9) Costamantc 2007; Reimcr ct al. 2008; (10) 
Anderhub et al. 2009b; Vittorini et al. 2009; (11) Abdo et al. 2011; (12) Abdo et al. 2010c; (13) Riiger et al. 2010; Abdo et al. 2009; (14) 
Tavecchio et al. 2010b; (15) Acciari et al. 2010c; (16) Acciari et al. 2010b; (17) Acciari et al. 2009c; Abdo et al. 2009; Tavecchio et al. 2010b; (18) 
Rugamcr ct al. 2011; Abdo ct al. 2009; (19) HESS Collaboration ct al. 2010b; (20) Aharonian ct al. 2007b; Abdo et al. 2010b; (21) Tavecchio et 
al. 2010b; (21) Aharonian et al. 2008; (23) Wolter et al. 2008 ; Abdo et al. 2009; (24) Aharonian et al. 2010. 
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logv [Hz] logv [Hz] 

Fig. 1. — Observed SEDs {scattered data points) with our best model fits (lines). Solid 
squares are for tlie data observed simultaneously and opened squares are for the data obtained 
not simultaneously. The data of high and low states are marked with blue and red symbols, 
respectively. If only one SED is obtained, the data are shown with black symbols. The 
Fermi/LAT observations are presented with a bow-ties or a magenta triangle (upper-limits 
for six sources, namely, lES 1101-232, lES 0347-121, H 2356-309, RGB J0152+017, lES 
0229-1-200, and PKS 0548-322). Dashed line of bow-tie indicates that the data are not used 
to calculate the values of with fitting results. 
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Fig. 1 — continued 
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Fig. 2. — Contours of the probability p in the B — S plane {panels a, d, g) along with the 
one-dimensional distributions of p for 5 and B {panels b, c, e, f, h, i). Gaussian or one-sided 
Gaussian function fits to the one-dimensional distributions are shown with thick solid lines 
{panels b, c, e, f, h, i). The values of p that correspond to the 1 a region of the Gaussian 
fits {panels h, c, e, and f) or the peak of the one-sided Gaussian function {panels h and i) 
are shown with dotted lines. The contours of la are marked with the thick lines in panels 
a and d, and the filled contour in Panel g is for p that corresponds to the peak value of the 
one-sided Gaussian function. The stars stand for the best fitting parameter set of S and B. 
The vertical dashed lines in panels a, d, g mark the lower limit of 5 derived from the 7-ray 
transparency condition. The dotted line in Panel a is the best fitting line for the 1 a region 
and the slope is ~ —2.33. 
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Fig. 3. — Distributions of the magnetic field B {Panel a), the beaming factor 6 {Panel b), the 
break Lorentz factor 7b of the electron distribution {Panel c), and the size of the radiating 
region {Panel d) for the 24 SEDs. Panel e — Distributions of the synchrotron radiation peak 
frequency z/g {red dashed line) and the inverse Compton scattering peak frequency Uc {blue 
solid line) for the 34 SEDs derived from our theoretical model fits. 
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[erg/s] log M^^M^ [erg/s] 

Fig. 4. — Panel a, b — Jet power as a function of the bolometric luminosity and the BH mass. 
Blue and red squares are for the sources in the high and low states, respectively, and black 
stars are for the sources with only one SED available. The up-triangle indicates that it is a 
lower limit of Pjet and no error of Lboi is available. The lines are the best fit lines (including 
the data of limits), which are log Pjet = (55.2 ± 3.1) — (1.21 ± 0.35) logMen for the low state 
data {dashed line) and log Pjet = (56.2 ± 4.1) - (1.22 ± 0.47) logMen for the high state data 
{solid line). Panel c — The powers associated with relativistic electrons (Pe) as a function of 
the powers of Poynting flux (Pg) in the jets. The line is the equality line. 
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Fig. 5. — Distributions of the powers associated with relativistic electrons Pe, cold protons 
Pp, Poynting flux Pg, radiation component P^, and the total power P^^t of the jets {left 
panel), together with the distributions of the ratios of these powers to the total jet power 
{right panels) for the 24 SEDs. 
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Fig. 6. — Comparison of the peak frequencies (and i/c) between the high and low states. 
Circles are for and triangles are for Uc. Since is a input parameter, no error of i/g is 
available. The solid line is the equality line. 
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Fig. 7. — Ratios of i?x = x^/x^ as a function of the ratio Ri TeV = -^rTev/-^iTeV> where 
x=[(5, B, Lboi, -Pjet, 7b], F is the flux density at 1 TeV, and the letters "H" and "L" mark 
the high and low states, respectively. The up-triangles and down-triangles indicate that 
they are the lower limits and the upper limits, and the opened circles indicate that no 
errors are available. The best fit line (including the data of limits) for R^^^ — Ri TeV is 
logi?^,, = (-0.13 ± 0.3) + (0.54 ± 0.27) logi?i TeV- 
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U' [erg cm"'] [Hz] 

Fig. 8. — Panel a — 7b as a function of energy density U' . Panel b and Panel c — Bolometric 
luminosity (Lboi) and the ratio of the two bump luminosities Ljc/i^syn as a function of 
synchrotron radiation peak frequency Ug. The gray circles in Panel a are the sample data 
from Ghisellini et al. (2010), other symbols with other color are the same as in Fig. HI and 
the diamonds indicate that those data points should move toward the left top of Figure. 
The best fit lines are logLic/^vgyn = (3.27± 1.0) — (0.23 ± 0.06) logz/g for the low state data 
{dashed line) and log Ljc /-^^syn = (3.5 ± 0.95) — (0.2 ± 0.06) log z/g for the high state data 
(excluding the source PKS 2005-489 in the high state; solid line). 
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